SD7213 
Keenan 

CLAIMS 

We claim: 

1 . A method for analyzing multivariate images, comprising: 

a) providing a data matrix D containing measured spectral data, 

b) transforming the data matrix D, using a wavelet transform, to obtain 

a transformed data matrix D , 

c) performing an Image analysis on the transformed data matrix D to 
obtain a transformed concentration matrix C and a spectral shapes matrix S, 
and 

d) computing a concentration matrix C from the transformed 
concentration matrix C . 

2. The method of Claim 1 , wherein the data matrix D comprises a total of j 
blocks of data D/, each data block D/ thereby providing a concentration block 0/ in 
step a), and wherein steps a) through d) are repeated sequentially until the 
concentration matrix C is accumulated blockwise, according to 

c = [c, c, ... c,., cj. 

3. The method of Claim 1 , wherein the wavelet transform comprises a Haar 
transform. 

4. The method of Claim 1 , further comprising thresholding the wavelet 
coefficients of the transformed data matrix D . 

5. The method of Claim 4, wherein the thresholding comprises decimating 
the detail coefficients. 

6. The method of Claim' 1 , wherein the image analysis of step c) comprises 
an alternating least squares analysis and the transformed concentration matrix 

C and the spectral shapes matrix S are obtained from a constrained least 
squares solution of min 5 -CS^ . 
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7. The method of Claim 6, wherein the altemating least squares analysis 
comprises a transformed non-negativity constraint. 

8. The method of Claim 1 , wherein the.computing step d) comprises applying 
an inverse wavelet transform to the transformed concentration matrix C to 
provide the concentration matrix C. 

9. The method of Claim 1 , wherein the computing step d) comprises 
projecting the data matrix D from step a) onto the spectral shapes matrix S from 

step c), according to minlb - CS^ . 



10. A method for analyzing multivariate images, comprising: 

a) providing a data factor matrix A and a data factor matrix B obtained 
from a factorization of measured spectral data D, 

b) transforming the data factor matrix A, using a wavelet transform, to 

obtain a transformed data factor matrix A , 

c) performing an image analysis on the transformed data factor matrix 

A and data factor matrix B to obtain a transformed concentration matrix C and 

a spectral shapes matrix S, and 

d) computing a concentration matrix 0 from the transformed 

concentration matrix C . 

1 1 . The method of Claim 1 0, wherein the data factor matrix A comprises a 
total of I blocks of data factors A,- and the data factor matrix B comprises k blocks 
of data factors B/, thereby providing a concentration block C, in step d), and 
wherein steps a) through d) are repeated sequentially until the concentration 
matrix C is accumulated blockwise, according to C = [c^ Cy.^ Cy] . 

12. The method of Claim 10, wherein the wavelet transform comprises a Haar 
transform. 

1 3. The method of Claim 1 0, further comprising thresholding the wavelet 
coefficients of the transformed data factor matrix A . 
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14. The method of Claim 13, wherein the thresholding comprises decimating 
the detail coefficients. 

15. The method of Claim 10, wherein the image analysis of step c) comprises 
an alternating least squares analysis and the transformed concentration matrix 

C and the spectral shapes matrix S are obtained from a constrained least 
squares solution of min AB^ -CS^ 

16. The method of Claim 15, wherein the alternating least squares analysis 
comprises a transformed non-negativity constraint. 

17. The method of Claim 10, wherein the computing step d) comprises 
applying an inverse wavelet transform to the transformed concentration matrix C 
to provide the concentration matrix C. 

18. The method of Claim 10, wherein the computing step d) comprises 
projecting the product of the data factor matrix A and the data factor matrix B 
from step a) onto the spectral shapes matrix S from step c), according to 

minllAB^ -CS^l and subject to appropriate constraints. 



19. The method of Claim 10, wherein the data factor matrix A comprises a 
scores matrix T and the data factor matrix B comprises a loadings matrix P, and 
wherein T and P are obtained from a principal components analysis of the 
measured spectral data D, according to D = TP^. 

20, The method of Claim 19, wherein T and P represent the significant 
components of the principal components. 
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